3.1.53 \(\int \frac {(c+d x)^2}{a+b \cot (e+f x)} \, dx\) [53]

3.1.53.1 Optimal result
3.1.53.2 Mathematica [A] (verified)
3.1.53.3 Rubi [A] (verified)
3.1.53.4 Maple [B] (verified)
3.1.53.5 Fricas [B] (verification not implemented)
3.1.53.6 Sympy [F]
3.1.53.7 Maxima [B] (verification not implemented)
3.1.53.8 Giac [F]
3.1.53.9 Mupad [F(-1)]

3.1.53.1 Optimal result

Integrand size = 20, antiderivative size = 181 \[ \int \frac {(c+d x)^2}{a+b \cot (e+f x)} \, dx=\frac {(c+d x)^3}{3 (a-i b) d}-\frac {b (c+d x)^2 \log \left (1-\frac {(a+i b) e^{2 i (e+f x)}}{a-i b}\right )}{\left (a^2+b^2\right ) f}+\frac {i b d (c+d x) \operatorname {PolyLog}\left (2,\frac {(a+i b) e^{2 i (e+f x)}}{a-i b}\right )}{\left (a^2+b^2\right ) f^2}-\frac {b d^2 \operatorname {PolyLog}\left (3,\frac {(a+i b) e^{2 i (e+f x)}}{a-i b}\right )}{2 \left (a^2+b^2\right ) f^3} \]

output
1/3*(d*x+c)^3/(a-I*b)/d-b*(d*x+c)^2*ln(1-(a+I*b)*exp(2*I*(f*x+e))/(a-I*b)) 
/(a^2+b^2)/f+I*b*d*(d*x+c)*polylog(2,(a+I*b)*exp(2*I*(f*x+e))/(a-I*b))/(a^ 
2+b^2)/f^2-1/2*b*d^2*polylog(3,(a+I*b)*exp(2*I*(f*x+e))/(a-I*b))/(a^2+b^2) 
/f^3
 
3.1.53.2 Mathematica [A] (verified)

Time = 1.64 (sec) , antiderivative size = 289, normalized size of antiderivative = 1.60 \[ \int \frac {(c+d x)^2}{a+b \cot (e+f x)} \, dx=\frac {b \left (\frac {4 i (c+d x)^3}{(a+i b) d}-\frac {6 \left (a \left (-1+e^{2 i e}\right )+i b \left (1+e^{2 i e}\right )\right ) (c+d x)^2 \log \left (1+\frac {(-a+i b) e^{-2 i (e+f x)}}{a+i b}\right )}{\left (a^2+b^2\right ) f}+\frac {3 d \left (-i a \left (-1+e^{2 i e}\right )+b \left (1+e^{2 i e}\right )\right ) \left (2 f (c+d x) \operatorname {PolyLog}\left (2,\frac {(a-i b) e^{-2 i (e+f x)}}{a+i b}\right )-i d \operatorname {PolyLog}\left (3,\frac {(a-i b) e^{-2 i (e+f x)}}{a+i b}\right )\right )}{\left (a^2+b^2\right ) f^3}\right )}{6 \left (a \left (-1+e^{2 i e}\right )+i b \left (1+e^{2 i e}\right )\right )}+\frac {x \left (3 c^2+3 c d x+d^2 x^2\right ) \sin (e)}{3 (b \cos (e)+a \sin (e))} \]

input
Integrate[(c + d*x)^2/(a + b*Cot[e + f*x]),x]
 
output
(b*(((4*I)*(c + d*x)^3)/((a + I*b)*d) - (6*(a*(-1 + E^((2*I)*e)) + I*b*(1 
+ E^((2*I)*e)))*(c + d*x)^2*Log[1 + (-a + I*b)/((a + I*b)*E^((2*I)*(e + f* 
x)))])/((a^2 + b^2)*f) + (3*d*((-I)*a*(-1 + E^((2*I)*e)) + b*(1 + E^((2*I) 
*e)))*(2*f*(c + d*x)*PolyLog[2, (a - I*b)/((a + I*b)*E^((2*I)*(e + f*x)))] 
 - I*d*PolyLog[3, (a - I*b)/((a + I*b)*E^((2*I)*(e + f*x)))]))/((a^2 + b^2 
)*f^3)))/(6*(a*(-1 + E^((2*I)*e)) + I*b*(1 + E^((2*I)*e)))) + (x*(3*c^2 + 
3*c*d*x + d^2*x^2)*Sin[e])/(3*(b*Cos[e] + a*Sin[e]))
 
3.1.53.3 Rubi [A] (verified)

Time = 0.70 (sec) , antiderivative size = 187, normalized size of antiderivative = 1.03, number of steps used = 8, number of rules used = 7, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.350, Rules used = {3042, 4214, 25, 2620, 3011, 2720, 7143}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int \frac {(c+d x)^2}{a+b \cot (e+f x)} \, dx\)

\(\Big \downarrow \) 3042

\(\displaystyle \int \frac {(c+d x)^2}{a-b \tan \left (e+f x+\frac {\pi }{2}\right )}dx\)

\(\Big \downarrow \) 4214

\(\displaystyle \frac {(c+d x)^3}{3 d (a-i b)}-2 i b \int -\frac {e^{2 i (e+f x)} (c+d x)^2}{(a-i b)^2-\left (a^2+b^2\right ) e^{2 i (e+f x)}}dx\)

\(\Big \downarrow \) 25

\(\displaystyle 2 i b \int \frac {e^{2 i (e+f x)} (c+d x)^2}{(a-i b)^2-\left (a^2+b^2\right ) e^{2 i (e+f x)}}dx+\frac {(c+d x)^3}{3 d (a-i b)}\)

\(\Big \downarrow \) 2620

\(\displaystyle 2 i b \left (\frac {i (c+d x)^2 \log \left (1-\frac {(a+i b) e^{2 i (e+f x)}}{a-i b}\right )}{2 f \left (a^2+b^2\right )}-\frac {i d \int (c+d x) \log \left (1-\frac {(a+i b) e^{2 i (e+f x)}}{a-i b}\right )dx}{f \left (a^2+b^2\right )}\right )+\frac {(c+d x)^3}{3 d (a-i b)}\)

\(\Big \downarrow \) 3011

\(\displaystyle 2 i b \left (\frac {i (c+d x)^2 \log \left (1-\frac {(a+i b) e^{2 i (e+f x)}}{a-i b}\right )}{2 f \left (a^2+b^2\right )}-\frac {i d \left (\frac {i (c+d x) \operatorname {PolyLog}\left (2,\frac {(a+i b) e^{2 i (e+f x)}}{a-i b}\right )}{2 f}-\frac {i d \int \operatorname {PolyLog}\left (2,\frac {(a+i b) e^{2 i (e+f x)}}{a-i b}\right )dx}{2 f}\right )}{f \left (a^2+b^2\right )}\right )+\frac {(c+d x)^3}{3 d (a-i b)}\)

\(\Big \downarrow \) 2720

\(\displaystyle 2 i b \left (\frac {i (c+d x)^2 \log \left (1-\frac {(a+i b) e^{2 i (e+f x)}}{a-i b}\right )}{2 f \left (a^2+b^2\right )}-\frac {i d \left (\frac {i (c+d x) \operatorname {PolyLog}\left (2,\frac {(a+i b) e^{2 i (e+f x)}}{a-i b}\right )}{2 f}-\frac {d \int e^{-2 i (e+f x)} \operatorname {PolyLog}\left (2,\frac {(a+i b) e^{2 i (e+f x)}}{a-i b}\right )de^{2 i (e+f x)}}{4 f^2}\right )}{f \left (a^2+b^2\right )}\right )+\frac {(c+d x)^3}{3 d (a-i b)}\)

\(\Big \downarrow \) 7143

\(\displaystyle 2 i b \left (\frac {i (c+d x)^2 \log \left (1-\frac {(a+i b) e^{2 i (e+f x)}}{a-i b}\right )}{2 f \left (a^2+b^2\right )}-\frac {i d \left (\frac {i (c+d x) \operatorname {PolyLog}\left (2,\frac {(a+i b) e^{2 i (e+f x)}}{a-i b}\right )}{2 f}-\frac {d \operatorname {PolyLog}\left (3,\frac {(a+i b) e^{2 i (e+f x)}}{a-i b}\right )}{4 f^2}\right )}{f \left (a^2+b^2\right )}\right )+\frac {(c+d x)^3}{3 d (a-i b)}\)

input
Int[(c + d*x)^2/(a + b*Cot[e + f*x]),x]
 
output
(c + d*x)^3/(3*(a - I*b)*d) + (2*I)*b*(((I/2)*(c + d*x)^2*Log[1 - ((a + I* 
b)*E^((2*I)*(e + f*x)))/(a - I*b)])/((a^2 + b^2)*f) - (I*d*(((I/2)*(c + d* 
x)*PolyLog[2, ((a + I*b)*E^((2*I)*(e + f*x)))/(a - I*b)])/f - (d*PolyLog[3 
, ((a + I*b)*E^((2*I)*(e + f*x)))/(a - I*b)])/(4*f^2)))/((a^2 + b^2)*f))
 

3.1.53.3.1 Defintions of rubi rules used

rule 25
Int[-(Fx_), x_Symbol] :> Simp[Identity[-1]   Int[Fx, x], x]
 

rule 2620
Int[(((F_)^((g_.)*((e_.) + (f_.)*(x_))))^(n_.)*((c_.) + (d_.)*(x_))^(m_.))/ 
((a_) + (b_.)*((F_)^((g_.)*((e_.) + (f_.)*(x_))))^(n_.)), x_Symbol] :> Simp 
[((c + d*x)^m/(b*f*g*n*Log[F]))*Log[1 + b*((F^(g*(e + f*x)))^n/a)], x] - Si 
mp[d*(m/(b*f*g*n*Log[F]))   Int[(c + d*x)^(m - 1)*Log[1 + b*((F^(g*(e + f*x 
)))^n/a)], x], x] /; FreeQ[{F, a, b, c, d, e, f, g, n}, x] && IGtQ[m, 0]
 

rule 2720
Int[u_, x_Symbol] :> With[{v = FunctionOfExponential[u, x]}, Simp[v/D[v, x] 
   Subst[Int[FunctionOfExponentialFunction[u, x]/x, x], x, v], x]] /; Funct 
ionOfExponentialQ[u, x] &&  !MatchQ[u, (w_)*((a_.)*(v_)^(n_))^(m_) /; FreeQ 
[{a, m, n}, x] && IntegerQ[m*n]] &&  !MatchQ[u, E^((c_.)*((a_.) + (b_.)*x)) 
*(F_)[v_] /; FreeQ[{a, b, c}, x] && InverseFunctionQ[F[x]]]
 

rule 3011
Int[Log[1 + (e_.)*((F_)^((c_.)*((a_.) + (b_.)*(x_))))^(n_.)]*((f_.) + (g_.) 
*(x_))^(m_.), x_Symbol] :> Simp[(-(f + g*x)^m)*(PolyLog[2, (-e)*(F^(c*(a + 
b*x)))^n]/(b*c*n*Log[F])), x] + Simp[g*(m/(b*c*n*Log[F]))   Int[(f + g*x)^( 
m - 1)*PolyLog[2, (-e)*(F^(c*(a + b*x)))^n], x], x] /; FreeQ[{F, a, b, c, e 
, f, g, n}, x] && GtQ[m, 0]
 

rule 3042
Int[u_, x_Symbol] :> Int[DeactivateTrig[u, x], x] /; FunctionOfTrigOfLinear 
Q[u, x]
 

rule 4214
Int[((c_.) + (d_.)*(x_))^(m_.)/((a_) + (b_.)*tan[(e_.) + Pi*(k_.) + (f_.)*( 
x_)]), x_Symbol] :> Simp[(c + d*x)^(m + 1)/(d*(m + 1)*(a + I*b)), x] + Simp 
[2*I*b   Int[(c + d*x)^m*E^(2*I*k*Pi)*(E^Simp[2*I*(e + f*x), x]/((a + I*b)^ 
2 + (a^2 + b^2)*E^(2*I*k*Pi)*E^Simp[2*I*(e + f*x), x])), x], x] /; FreeQ[{a 
, b, c, d, e, f}, x] && IntegerQ[4*k] && NeQ[a^2 + b^2, 0] && IGtQ[m, 0]
 

rule 7143
Int[PolyLog[n_, (c_.)*((a_.) + (b_.)*(x_))^(p_.)]/((d_.) + (e_.)*(x_)), x_S 
ymbol] :> Simp[PolyLog[n + 1, c*(a + b*x)^p]/(e*p), x] /; FreeQ[{a, b, c, d 
, e, n, p}, x] && EqQ[b*d, a*e]
 
3.1.53.4 Maple [B] (verified)

Both result and optimal contain complex but leaf count of result is larger than twice the leaf count of optimal. 896 vs. \(2 (163 ) = 326\).

Time = 0.32 (sec) , antiderivative size = 897, normalized size of antiderivative = 4.96

method result size
risch \(\frac {d^{2} x^{3}}{3 i b +3 a}+\frac {d c \,x^{2}}{i b +a}+\frac {c^{2} x}{i b +a}+\frac {c^{3}}{3 d \left (i b +a \right )}+\frac {2 b \,d^{2} x^{3}}{3 \left (-i a +b \right ) \left (-i b +a \right )}+\frac {2 i b \,c^{2} \ln \left ({\mathrm e}^{i \left (f x +e \right )}\right )}{f \left (-i a +b \right ) \left (i b -a \right )}-\frac {i b \,e^{2} d^{2} \ln \left (1-\frac {\left (i b +a \right ) {\mathrm e}^{2 i \left (f x +e \right )}}{-i b +a}\right )}{f^{3} \left (-i a +b \right ) \left (-i b +a \right )}+\frac {2 i b c d \ln \left (1-\frac {\left (i b +a \right ) {\mathrm e}^{2 i \left (f x +e \right )}}{-i b +a}\right ) e}{f^{2} \left (-i a +b \right ) \left (-i b +a \right )}-\frac {i b \,e^{2} d^{2} \ln \left ({\mathrm e}^{2 i \left (f x +e \right )} a +i b \,{\mathrm e}^{2 i \left (f x +e \right )}-a +i b \right )}{f^{3} \left (-i a +b \right ) \left (i b -a \right )}-\frac {4 i b e c d \ln \left ({\mathrm e}^{i \left (f x +e \right )}\right )}{f^{2} \left (-i a +b \right ) \left (i b -a \right )}-\frac {2 b \,e^{2} d^{2} x}{f^{2} \left (-i a +b \right ) \left (-i b +a \right )}+\frac {b \,d^{2} \operatorname {polylog}\left (2, \frac {\left (i b +a \right ) {\mathrm e}^{2 i \left (f x +e \right )}}{-i b +a}\right ) x}{f^{2} \left (-i a +b \right ) \left (-i b +a \right )}+\frac {2 i b \,e^{2} d^{2} \ln \left ({\mathrm e}^{i \left (f x +e \right )}\right )}{f^{3} \left (-i a +b \right ) \left (i b -a \right )}+\frac {2 i b e c d \ln \left ({\mathrm e}^{2 i \left (f x +e \right )} a +i b \,{\mathrm e}^{2 i \left (f x +e \right )}-a +i b \right )}{f^{2} \left (-i a +b \right ) \left (i b -a \right )}+\frac {2 b c d \,x^{2}}{\left (-i a +b \right ) \left (-i b +a \right )}+\frac {4 b c d e x}{f \left (-i a +b \right ) \left (-i b +a \right )}+\frac {2 b c d \,e^{2}}{f^{2} \left (-i a +b \right ) \left (-i b +a \right )}+\frac {b c d \operatorname {polylog}\left (2, \frac {\left (i b +a \right ) {\mathrm e}^{2 i \left (f x +e \right )}}{-i b +a}\right )}{f^{2} \left (-i a +b \right ) \left (-i b +a \right )}+\frac {i b \,d^{2} \operatorname {polylog}\left (3, \frac {\left (i b +a \right ) {\mathrm e}^{2 i \left (f x +e \right )}}{-i b +a}\right )}{2 f^{3} \left (-i a +b \right ) \left (-i b +a \right )}-\frac {i b \,c^{2} \ln \left ({\mathrm e}^{2 i \left (f x +e \right )} a +i b \,{\mathrm e}^{2 i \left (f x +e \right )}-a +i b \right )}{f \left (-i a +b \right ) \left (i b -a \right )}+\frac {2 i b c d \ln \left (1-\frac {\left (i b +a \right ) {\mathrm e}^{2 i \left (f x +e \right )}}{-i b +a}\right ) x}{f \left (-i a +b \right ) \left (-i b +a \right )}-\frac {4 b \,e^{3} d^{2}}{3 f^{3} \left (-i a +b \right ) \left (-i b +a \right )}+\frac {i b \,d^{2} \ln \left (1-\frac {\left (i b +a \right ) {\mathrm e}^{2 i \left (f x +e \right )}}{-i b +a}\right ) x^{2}}{f \left (-i a +b \right ) \left (-i b +a \right )}\) \(897\)

input
int((d*x+c)^2/(a+b*cot(f*x+e)),x,method=_RETURNVERBOSE)
 
output
1/3*d^2/(a+I*b)*x^3+d/(a+I*b)*c*x^2+1/(a+I*b)*c^2*x+1/3/d/(a+I*b)*c^3+2/3/ 
(b-I*a)*b/(a-I*b)*d^2*x^3+2*I/f/(b-I*a)*b*c^2/(I*b-a)*ln(exp(I*(f*x+e)))-I 
/f^3/(b-I*a)*b/(a-I*b)*e^2*d^2*ln(1-(a+I*b)*exp(2*I*(f*x+e))/(a-I*b))+2*I/ 
f^2/(b-I*a)*b/(a-I*b)*c*d*ln(1-(a+I*b)*exp(2*I*(f*x+e))/(a-I*b))*e-I/f^3/( 
b-I*a)*b*e^2*d^2/(I*b-a)*ln(exp(2*I*(f*x+e))*a+I*b*exp(2*I*(f*x+e))-a+I*b) 
-4*I/f^2/(b-I*a)*b*e*c*d/(I*b-a)*ln(exp(I*(f*x+e)))-2/f^2/(b-I*a)*b/(a-I*b 
)*e^2*d^2*x+1/f^2/(b-I*a)*b/(a-I*b)*d^2*polylog(2,(a+I*b)*exp(2*I*(f*x+e)) 
/(a-I*b))*x+2*I/f^3/(b-I*a)*b*e^2*d^2/(I*b-a)*ln(exp(I*(f*x+e)))+2*I/f^2/( 
b-I*a)*b*e*c*d/(I*b-a)*ln(exp(2*I*(f*x+e))*a+I*b*exp(2*I*(f*x+e))-a+I*b)+2 
/(b-I*a)*b/(a-I*b)*c*d*x^2+4/f/(b-I*a)*b/(a-I*b)*c*d*e*x+2/f^2/(b-I*a)*b/( 
a-I*b)*c*d*e^2+1/f^2/(b-I*a)*b/(a-I*b)*c*d*polylog(2,(a+I*b)*exp(2*I*(f*x+ 
e))/(a-I*b))+1/2*I/f^3/(b-I*a)*b/(a-I*b)*d^2*polylog(3,(a+I*b)*exp(2*I*(f* 
x+e))/(a-I*b))-I/f/(b-I*a)*b*c^2/(I*b-a)*ln(exp(2*I*(f*x+e))*a+I*b*exp(2*I 
*(f*x+e))-a+I*b)+2*I/f/(b-I*a)*b/(a-I*b)*c*d*ln(1-(a+I*b)*exp(2*I*(f*x+e)) 
/(a-I*b))*x-4/3/f^3/(b-I*a)*b/(a-I*b)*e^3*d^2+I/f/(b-I*a)*b/(a-I*b)*d^2*ln 
(1-(a+I*b)*exp(2*I*(f*x+e))/(a-I*b))*x^2
 
3.1.53.5 Fricas [B] (verification not implemented)

Both result and optimal contain complex but leaf count of result is larger than twice the leaf count of optimal. 732 vs. \(2 (154) = 308\).

Time = 0.31 (sec) , antiderivative size = 732, normalized size of antiderivative = 4.04 \[ \int \frac {(c+d x)^2}{a+b \cot (e+f x)} \, dx=\frac {4 \, a d^{2} f^{3} x^{3} + 12 \, a c d f^{3} x^{2} + 12 \, a c^{2} f^{3} x - 3 \, b d^{2} {\rm polylog}\left (3, \frac {{\left (a^{2} + 2 i \, a b - b^{2}\right )} \cos \left (2 \, f x + 2 \, e\right ) + {\left (i \, a^{2} - 2 \, a b - i \, b^{2}\right )} \sin \left (2 \, f x + 2 \, e\right )}{a^{2} + b^{2}}\right ) - 3 \, b d^{2} {\rm polylog}\left (3, \frac {{\left (a^{2} - 2 i \, a b - b^{2}\right )} \cos \left (2 \, f x + 2 \, e\right ) + {\left (-i \, a^{2} - 2 \, a b + i \, b^{2}\right )} \sin \left (2 \, f x + 2 \, e\right )}{a^{2} + b^{2}}\right ) - 6 \, {\left (-i \, b d^{2} f x - i \, b c d f\right )} {\rm Li}_2\left (-\frac {a^{2} + b^{2} - {\left (a^{2} + 2 i \, a b - b^{2}\right )} \cos \left (2 \, f x + 2 \, e\right ) + {\left (-i \, a^{2} + 2 \, a b + i \, b^{2}\right )} \sin \left (2 \, f x + 2 \, e\right )}{a^{2} + b^{2}} + 1\right ) - 6 \, {\left (i \, b d^{2} f x + i \, b c d f\right )} {\rm Li}_2\left (-\frac {a^{2} + b^{2} - {\left (a^{2} - 2 i \, a b - b^{2}\right )} \cos \left (2 \, f x + 2 \, e\right ) + {\left (i \, a^{2} + 2 \, a b - i \, b^{2}\right )} \sin \left (2 \, f x + 2 \, e\right )}{a^{2} + b^{2}} + 1\right ) - 6 \, {\left (b d^{2} e^{2} - 2 \, b c d e f + b c^{2} f^{2}\right )} \log \left (\frac {1}{2} \, a^{2} + i \, a b - \frac {1}{2} \, b^{2} - \frac {1}{2} \, {\left (a^{2} + b^{2}\right )} \cos \left (2 \, f x + 2 \, e\right ) + \frac {1}{2} \, {\left (i \, a^{2} + i \, b^{2}\right )} \sin \left (2 \, f x + 2 \, e\right )\right ) - 6 \, {\left (b d^{2} e^{2} - 2 \, b c d e f + b c^{2} f^{2}\right )} \log \left (-\frac {1}{2} \, a^{2} + i \, a b + \frac {1}{2} \, b^{2} + \frac {1}{2} \, {\left (a^{2} + b^{2}\right )} \cos \left (2 \, f x + 2 \, e\right ) + \frac {1}{2} \, {\left (i \, a^{2} + i \, b^{2}\right )} \sin \left (2 \, f x + 2 \, e\right )\right ) - 6 \, {\left (b d^{2} f^{2} x^{2} + 2 \, b c d f^{2} x - b d^{2} e^{2} + 2 \, b c d e f\right )} \log \left (\frac {a^{2} + b^{2} - {\left (a^{2} + 2 i \, a b - b^{2}\right )} \cos \left (2 \, f x + 2 \, e\right ) + {\left (-i \, a^{2} + 2 \, a b + i \, b^{2}\right )} \sin \left (2 \, f x + 2 \, e\right )}{a^{2} + b^{2}}\right ) - 6 \, {\left (b d^{2} f^{2} x^{2} + 2 \, b c d f^{2} x - b d^{2} e^{2} + 2 \, b c d e f\right )} \log \left (\frac {a^{2} + b^{2} - {\left (a^{2} - 2 i \, a b - b^{2}\right )} \cos \left (2 \, f x + 2 \, e\right ) + {\left (i \, a^{2} + 2 \, a b - i \, b^{2}\right )} \sin \left (2 \, f x + 2 \, e\right )}{a^{2} + b^{2}}\right )}{12 \, {\left (a^{2} + b^{2}\right )} f^{3}} \]

input
integrate((d*x+c)^2/(a+b*cot(f*x+e)),x, algorithm="fricas")
 
output
1/12*(4*a*d^2*f^3*x^3 + 12*a*c*d*f^3*x^2 + 12*a*c^2*f^3*x - 3*b*d^2*polylo 
g(3, ((a^2 + 2*I*a*b - b^2)*cos(2*f*x + 2*e) + (I*a^2 - 2*a*b - I*b^2)*sin 
(2*f*x + 2*e))/(a^2 + b^2)) - 3*b*d^2*polylog(3, ((a^2 - 2*I*a*b - b^2)*co 
s(2*f*x + 2*e) + (-I*a^2 - 2*a*b + I*b^2)*sin(2*f*x + 2*e))/(a^2 + b^2)) - 
 6*(-I*b*d^2*f*x - I*b*c*d*f)*dilog(-(a^2 + b^2 - (a^2 + 2*I*a*b - b^2)*co 
s(2*f*x + 2*e) + (-I*a^2 + 2*a*b + I*b^2)*sin(2*f*x + 2*e))/(a^2 + b^2) + 
1) - 6*(I*b*d^2*f*x + I*b*c*d*f)*dilog(-(a^2 + b^2 - (a^2 - 2*I*a*b - b^2) 
*cos(2*f*x + 2*e) + (I*a^2 + 2*a*b - I*b^2)*sin(2*f*x + 2*e))/(a^2 + b^2) 
+ 1) - 6*(b*d^2*e^2 - 2*b*c*d*e*f + b*c^2*f^2)*log(1/2*a^2 + I*a*b - 1/2*b 
^2 - 1/2*(a^2 + b^2)*cos(2*f*x + 2*e) + 1/2*(I*a^2 + I*b^2)*sin(2*f*x + 2* 
e)) - 6*(b*d^2*e^2 - 2*b*c*d*e*f + b*c^2*f^2)*log(-1/2*a^2 + I*a*b + 1/2*b 
^2 + 1/2*(a^2 + b^2)*cos(2*f*x + 2*e) + 1/2*(I*a^2 + I*b^2)*sin(2*f*x + 2* 
e)) - 6*(b*d^2*f^2*x^2 + 2*b*c*d*f^2*x - b*d^2*e^2 + 2*b*c*d*e*f)*log((a^2 
 + b^2 - (a^2 + 2*I*a*b - b^2)*cos(2*f*x + 2*e) + (-I*a^2 + 2*a*b + I*b^2) 
*sin(2*f*x + 2*e))/(a^2 + b^2)) - 6*(b*d^2*f^2*x^2 + 2*b*c*d*f^2*x - b*d^2 
*e^2 + 2*b*c*d*e*f)*log((a^2 + b^2 - (a^2 - 2*I*a*b - b^2)*cos(2*f*x + 2*e 
) + (I*a^2 + 2*a*b - I*b^2)*sin(2*f*x + 2*e))/(a^2 + b^2)))/((a^2 + b^2)*f 
^3)
 
3.1.53.6 Sympy [F]

\[ \int \frac {(c+d x)^2}{a+b \cot (e+f x)} \, dx=\int \frac {\left (c + d x\right )^{2}}{a + b \cot {\left (e + f x \right )}}\, dx \]

input
integrate((d*x+c)**2/(a+b*cot(f*x+e)),x)
 
output
Integral((c + d*x)**2/(a + b*cot(e + f*x)), x)
 
3.1.53.7 Maxima [B] (verification not implemented)

Both result and optimal contain complex but leaf count of result is larger than twice the leaf count of optimal. 720 vs. \(2 (154) = 308\).

Time = 0.54 (sec) , antiderivative size = 720, normalized size of antiderivative = 3.98 \[ \int \frac {(c+d x)^2}{a+b \cot (e+f x)} \, dx=-\frac {6 \, c d e {\left (\frac {2 \, {\left (f x + e\right )} a}{{\left (a^{2} + b^{2}\right )} f} - \frac {2 \, b \log \left (a \tan \left (f x + e\right ) + b\right )}{{\left (a^{2} + b^{2}\right )} f} + \frac {b \log \left (\tan \left (f x + e\right )^{2} + 1\right )}{{\left (a^{2} + b^{2}\right )} f}\right )} - 3 \, {\left (\frac {2 \, {\left (f x + e\right )} a}{a^{2} + b^{2}} - \frac {2 \, b \log \left (a \tan \left (f x + e\right ) + b\right )}{a^{2} + b^{2}} + \frac {b \log \left (\tan \left (f x + e\right )^{2} + 1\right )}{a^{2} + b^{2}}\right )} c^{2} - \frac {2 \, {\left (f x + e\right )}^{3} {\left (a + i \, b\right )} d^{2} + 6 \, {\left (f x + e\right )} {\left (a + i \, b\right )} d^{2} e^{2} - 6 i \, b d^{2} e^{2} \arctan \left (b \cos \left (2 \, f x + 2 \, e\right ) + a \sin \left (2 \, f x + 2 \, e\right ) + b, a \cos \left (2 \, f x + 2 \, e\right ) - b \sin \left (2 \, f x + 2 \, e\right ) - a\right ) - 3 \, b d^{2} e^{2} \log \left ({\left (a^{2} + b^{2}\right )} \cos \left (2 \, f x + 2 \, e\right )^{2} + 4 \, a b \sin \left (2 \, f x + 2 \, e\right ) + {\left (a^{2} + b^{2}\right )} \sin \left (2 \, f x + 2 \, e\right )^{2} + a^{2} + b^{2} - 2 \, {\left (a^{2} - b^{2}\right )} \cos \left (2 \, f x + 2 \, e\right )\right ) - 3 \, b d^{2} {\rm Li}_{3}(\frac {{\left (i \, a - b\right )} e^{\left (2 i \, f x + 2 i \, e\right )}}{i \, a + b}) - 6 \, {\left ({\left (a + i \, b\right )} d^{2} e - {\left (a + i \, b\right )} c d f\right )} {\left (f x + e\right )}^{2} + 6 \, {\left (-i \, {\left (f x + e\right )}^{2} b d^{2} + 2 \, {\left (i \, b d^{2} e - i \, b c d f\right )} {\left (f x + e\right )}\right )} \arctan \left (-\frac {2 \, a b \cos \left (2 \, f x + 2 \, e\right ) + {\left (a^{2} - b^{2}\right )} \sin \left (2 \, f x + 2 \, e\right )}{a^{2} + b^{2}}, \frac {2 \, a b \sin \left (2 \, f x + 2 \, e\right ) + a^{2} + b^{2} - {\left (a^{2} - b^{2}\right )} \cos \left (2 \, f x + 2 \, e\right )}{a^{2} + b^{2}}\right ) + 6 \, {\left (i \, {\left (f x + e\right )} b d^{2} - i \, b d^{2} e + i \, b c d f\right )} {\rm Li}_2\left (\frac {{\left (i \, a - b\right )} e^{\left (2 i \, f x + 2 i \, e\right )}}{i \, a + b}\right ) - 3 \, {\left ({\left (f x + e\right )}^{2} b d^{2} - 2 \, {\left (b d^{2} e - b c d f\right )} {\left (f x + e\right )}\right )} \log \left (\frac {{\left (a^{2} + b^{2}\right )} \cos \left (2 \, f x + 2 \, e\right )^{2} + 4 \, a b \sin \left (2 \, f x + 2 \, e\right ) + {\left (a^{2} + b^{2}\right )} \sin \left (2 \, f x + 2 \, e\right )^{2} + a^{2} + b^{2} - 2 \, {\left (a^{2} - b^{2}\right )} \cos \left (2 \, f x + 2 \, e\right )}{a^{2} + b^{2}}\right )}{{\left (a^{2} + b^{2}\right )} f^{2}}}{6 \, f} \]

input
integrate((d*x+c)^2/(a+b*cot(f*x+e)),x, algorithm="maxima")
 
output
-1/6*(6*c*d*e*(2*(f*x + e)*a/((a^2 + b^2)*f) - 2*b*log(a*tan(f*x + e) + b) 
/((a^2 + b^2)*f) + b*log(tan(f*x + e)^2 + 1)/((a^2 + b^2)*f)) - 3*(2*(f*x 
+ e)*a/(a^2 + b^2) - 2*b*log(a*tan(f*x + e) + b)/(a^2 + b^2) + b*log(tan(f 
*x + e)^2 + 1)/(a^2 + b^2))*c^2 - (2*(f*x + e)^3*(a + I*b)*d^2 + 6*(f*x + 
e)*(a + I*b)*d^2*e^2 - 6*I*b*d^2*e^2*arctan2(b*cos(2*f*x + 2*e) + a*sin(2* 
f*x + 2*e) + b, a*cos(2*f*x + 2*e) - b*sin(2*f*x + 2*e) - a) - 3*b*d^2*e^2 
*log((a^2 + b^2)*cos(2*f*x + 2*e)^2 + 4*a*b*sin(2*f*x + 2*e) + (a^2 + b^2) 
*sin(2*f*x + 2*e)^2 + a^2 + b^2 - 2*(a^2 - b^2)*cos(2*f*x + 2*e)) - 3*b*d^ 
2*polylog(3, (I*a - b)*e^(2*I*f*x + 2*I*e)/(I*a + b)) - 6*((a + I*b)*d^2*e 
 - (a + I*b)*c*d*f)*(f*x + e)^2 + 6*(-I*(f*x + e)^2*b*d^2 + 2*(I*b*d^2*e - 
 I*b*c*d*f)*(f*x + e))*arctan2(-(2*a*b*cos(2*f*x + 2*e) + (a^2 - b^2)*sin( 
2*f*x + 2*e))/(a^2 + b^2), (2*a*b*sin(2*f*x + 2*e) + a^2 + b^2 - (a^2 - b^ 
2)*cos(2*f*x + 2*e))/(a^2 + b^2)) + 6*(I*(f*x + e)*b*d^2 - I*b*d^2*e + I*b 
*c*d*f)*dilog((I*a - b)*e^(2*I*f*x + 2*I*e)/(I*a + b)) - 3*((f*x + e)^2*b* 
d^2 - 2*(b*d^2*e - b*c*d*f)*(f*x + e))*log(((a^2 + b^2)*cos(2*f*x + 2*e)^2 
 + 4*a*b*sin(2*f*x + 2*e) + (a^2 + b^2)*sin(2*f*x + 2*e)^2 + a^2 + b^2 - 2 
*(a^2 - b^2)*cos(2*f*x + 2*e))/(a^2 + b^2)))/((a^2 + b^2)*f^2))/f
 
3.1.53.8 Giac [F]

\[ \int \frac {(c+d x)^2}{a+b \cot (e+f x)} \, dx=\int { \frac {{\left (d x + c\right )}^{2}}{b \cot \left (f x + e\right ) + a} \,d x } \]

input
integrate((d*x+c)^2/(a+b*cot(f*x+e)),x, algorithm="giac")
 
output
integrate((d*x + c)^2/(b*cot(f*x + e) + a), x)
 
3.1.53.9 Mupad [F(-1)]

Timed out. \[ \int \frac {(c+d x)^2}{a+b \cot (e+f x)} \, dx=\int \frac {{\left (c+d\,x\right )}^2}{a+b\,\mathrm {cot}\left (e+f\,x\right )} \,d x \]

input
int((c + d*x)^2/(a + b*cot(e + f*x)),x)
 
output
int((c + d*x)^2/(a + b*cot(e + f*x)), x)